Interplay between chiral and deconfinement phase transitions 
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Abstract. By using the dressed Polyakov loop or dual chiral condensate as an equivalent order parameter 
of the deconfinement phase transition, we investigate the relation between the chiral and deconfinement phase 
transitions at finite temperature and density in the framework of three-flavor Nambu-Jona-Lasinio (NJL) model. 
It is found that in the chiral limit, the critical temperature for chiral phase transition coincides with that of the 
dressed Polyakov loop in the whole (T,pi) plane. In the case of explicit chiral symmetry breaking, it is found that 
the phase transitions are flavor dependent. For each flavor, the transition temperature for chiral restoration T* 
is smaller than that of the dressed Polyakov loop in the low baryon density region where the transition is a 
crossover, and, the two critical temperatures coincide in the high baryon density region where the phase transition 
is of first order. Therefore, there are two critical end points, i.e, T^ Ep and T S CEP at finite density. We also explain 
the feature of T x c = T® in the case of 1st and 2nd order phase transitions, and T* < T® in the case of crossover, 
and expect this feature is general and can be extended to full QCD theory. 



1 Introduction 

The interplay between chiral symmetry breaking and con- 
finement as well as the chiral and deconfinement phase 
transitions at finite temperature and density are of contin- 
uous interests [1 2 3 4 5I6I7I8I9I . The two transitions are 
characterized by the breaking and restoration of chiral and 
center symmetry, which are well defined in two extreme 
quark mass limits, respectively. In the chiral limit when 
the current quark mass is zero m = 0, the chiral condensate 
is the order parameter for the chiral phase transition. When 
the current quark mass goes to infinity m — > oo, QCD be- 
comes pure gauge S U(3) theory, which is center symmet- 
ric in the vacuum, and the usually used order parameter is 
the Polyakov loop (TJ. 

The relation between the chiral and deconfinement phase 
transitions has attracted more interest recently in studying 
the phase diagram at high baryon density region IflOl . It 
is conjectured in Ref. ifTTI that in large N c limit, a con- 
fined but chiral symmetric phase, which is called quarky- 
onic phase can exist in the high baryon density region. It is 
very interesting to study whether this quarkyonic phase can 
survive in real QCD phase diagram, and how it competes 
with nuclear matter and the color superconducting phase 
02) • (However, it is worthy of noticing that in Ref. [13|, 
it is found that at zero chemical potential, the lattice re- 
sults for the thermodynamical properties have a very mild 
dependence on the number of colors.) 

Lattice QCD at the current stage cannot go to very 
high baryon density. For zero chemical potential, previ- 
ous lattice results show that the chiral and deconfinement 
phase transitions occur at the same temperature, e.g, in 
Ref. 1141151161171181 . and also in review papers [19 20 1. 



In recent years, three lattice groups, MILC group ETI . 
RBC-Bielefeld group 11221 and Wuppertal-Budapest group 
[23 24 25 1 have studied the chiral and deconfinement phase 
transition temperatures with almost physical quark masses. 
The RBC-Bielefeld group claimed that the two critical tem- 
peratures for Nf = 2 + 1 coincide at T c = 192(7)(4)MeV. 
The Wuppetal-Budapest group found that for the case of 
Nf = 2 + 1, there are three transition temperatures, the 
transition temperature for chiral restoration of u, d quarks 
jx(ud) _ 151 (3)(3) MeV, me transition temperature forchi- 
ral restoration of s quark r* (s) = 175(2)(4)MeV and the de- 
confinement transition temperature Tf. - 176(3)(4)MeV. 
From this result, we can read that the critical tempera- 
tures for different transitions are different. According to 
the Wuppetal-Budapest group, this is the consequence of 
the crossover nature. 

In the framework of QCD effective models, there is 
still no dynamical model which can describe the chiral 
symmetry breaking and confinement simultaneously. The 
main difficulty of effective QCD model to include confine- 
ment mechanism lies in that it is difficult to calculate the 
Polyakov loop analytically. Currently, the popular mod- 
els used to investigate the chiral and deconfinement phase 
transitions are the Polyakov Nambu-Jona-Lasinio model 
(PNJL) I26I27I28I29I30I3 11321331 and Polyakov linear sigma 
model (PLSM) 0341351 . However, the shortcoming of these 
models is that the temperature dependence of the Polyakov- 
loop potential is put in by hand from lattice result, which 
cannot be self-consistently extended to finite baryon den- 
sity. Recently, efforts have been made in Ref.[36| to derive 
a low-energy effective theory for confinement-deconfinement 
and chiral-symmetry breaking/restoration from first-principle. 
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Recent investigation revealed that quark propagator, heat 
kernels can also act as an order parameter as they transform 
non trivially under the center transformation related to de- 
confinement transition [37 38 39 1. But the exciting result 
is the behavior of spectral sum of the Dirac operator under 
center transformation 13 814014 11421 . A new order parame- 
ter, called dressed Polyakov loop has been defined which 
can be represented as a spectral sum of the Dirac operator 
|42|. It has been found the infrared part of the spectrum 
particularly plays a leading role in confinement [38]. This 
result is encouraging since it gives a hope to relate the chi- 
ral phase transition with the confinement-deconfinement 
phase transition. The order parameter for chiral phase tran- 
sition is related to the spectral density of the Dirac opera- 
tor through Banks-Casher relation [4 |. Therefore, both the 
dressed Polyakov loop and the chiral condensate are re- 
lated to the spectral sum of the Dirac operator. 

Behavior of the dressed Polyakov loop is mainly stud- 
ied in the framework of Lattice gauge theory 143 1441 . Apart 
from that, studies based on Dyson-Schwinger equations 
1451461471 and PNJL model [48 1 have been carried out. In 
those studies the role of dressed Polyakov loop as an order 
parameter is discussed at zero chemical potential. In this 
paper, we show our results 11491 of investigating the QCD 
phase diagram at finite temperature and density by using 
the dressed Polyakov loop as an equivalent order parame- 
ter in the Nambu-Jona-Lasinio (NJL) model . It is known 
that he NJL model lacks of confinement and the gluon dy- 
namics is encoded in a static coupling constant for four 
point contact interaction. However, assuming that we can 
read the information of confinement from the dual chiral 
condensate, it would be interesting to see the behavior of 
the dressed Polyakov loop in a scenario without any ex- 
plicit mechanism for confinement. 

In this paper, we show the phase transitions in the two- 
flavor and three-flavor NJL models in {T,jj) plane in chi- 
ral limit as well for small quark mass limit. This paper is 
organized as follows: We introduce the dressed Polyakov 
loop as an equivalent order parameter of confinement de- 
confinement phase transition and the NJL model in Sec. 
|2] Then in Secf3] we show the results of two-flavor QCD 
phase diagram in T — fi plane in the chiral limit and in 
the case of explicit chiral symmetry breaking, respectively. 
We offer an analysis on the relation between the chiral and 
deconfinement phase transitions in Sec. [4] In Sec|5] we 
show the phase diagram at finite temperature and density 
for three-flavor case. At the end, we give the conclusion 
and discussion. 



2 Dressed Polyakov loop and the NJL 
model 

We firstly introduce the dressed Polyakov loop. To do this 
we have to consider a U( 1 ) valued boundary condition for 
the fermionic fields in the temporal direction instead of the 
canonical choice of anti-periodic boundary condition, 



where < <f> < 2n is the phase angle and (3 is the inverse 
temperature. 

Dual quark condensate E n is then defined by the Fourier 
transform (w.r.t the phase </>) of the general boundary con- 
dition dependent quark condensate [42 43 44 1, 
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where n is the winding number. 

Particular case of n = 1 is called the dressed Polyakov 
loop which transforms in the same way as the conven- 
tional thin Polyakov loop under the center symmetry and 
hence is an order parameter for the deconfinement transi- 
tion [42 43 44 1 . It reduces to the thin Polyakov loop and 
to the dual of the conventional chiral condensate in infinite 
and zero quark mass limits respectively, i.e., in the chiral 
limit m — > we get the dual of the conventional chiral 
condensate and in the m — > oo limit we have thin Polyakov 
loop 1421431441 . 

The Lagrangian of three-flavor NJL model l50l is given 

as 

L = <A0y<9„ - m)^ + G s 2 {(<A^) 2 + (jpiysr^) 2 } 

a 

- ^Det/t^l + 75 )iff] + Det0(l - y 5 ¥]}. (3) 

Where if/ = (u, d, s) T denotes the transpose of the quark 
field, and m = Diag(m„, m^, m s ) is the corresponding mass 
matrix in the flavor space. T a with a = 1, • • • , Nj : - 1 are 
the eight Gell-Mann matrices, and Dety means determi- 
nant in flavor space. The last term is the standard form of 
the 't Hooft interaction, which is invariant under S U(3)l X 
SU{3)r X U(1)b symmetry, but breaks down the Ua(1) 
symmetry. 

The <p dependent thermodynamic potential in the mean 
field level is given as following: 

= ^ Q^M f + 2G S 2 <°">1/ - 4 ^<cr>0,«< cr >0,<<<cr>0i>4) 

/ / 



with 
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Where the sum is in the flavor space, E p j = Jp 2 + 
and E ± , = E p j + [p.+i(<f>-7f)T], with the constituent quark 



mass 



iff(x,P) = e-'VO,0), 



(1) 



= mi - AG s {cr)0,i + IKia-^jicr)^, (6) 
where (i, j, k) is the quark flavor indices ( u, d, s), and {cr)<pj - 
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3 Phase diagram for two-flavor case 

We firstly show the results in the two-flavor case and con- 
sider the isospin symmetric limit, i.e, we take Nf = 2, 
K = and m u = in Eq.©. The thermodynamic po- 
tential contains imaginary part. We take only the real part 
of the potential and the imaginary phase factor is not con- 
sidered in this work. The mean field (cr)^ is obtained by 
minimizing the potential for each value of e [0, 2n) 
for fixed T and /i. The conventional chiral condensate is 
(o")tt = (ifrifdn- The dressed Polyakov loop E\ is obtained 
by integrating over the angle. The values of the parame- 
ters A and G s are taken as 0.6315GeV and 5.498GeV~ 2 , 
respectively. 

We investigate phase transitions for two cases, i.e., in 
the chiral limit and with explicit chiral symmetry breaking. 
Fig- Hand [2] show the behavior of the angle dependence 
of the general chiral condensate for various chemical po- 
tentials and temperatures for m = and m - 5.5MeV, 
respectively. The four curves presented in each figure rep- 
resent two temperatures above and below the critical tem- 
perature for two particular values of the chemical poten- 
tial. Same qualitative features have been found for both the 
quark masses. The variation is symmetrical around (f> — n 
as reported in other studies H45I48H . 
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Fig. 2. Angle variation of (cr)^ for different values of tempera- 
tures and chemical potentials in the case of m = 5.5MeV. The 
solid line corresponds to T = 150MeV,/j = lOOMeV, dashed 
line corresponds to T = 250MeV,/j = lOOMeV, dash-dotted line 
corresponds to T = 20MeV, \i = 340MeV and dotted line corre- 
sponds to T = 40MeV iy u = 340MeV. 

is consistent with the expectation of complete restoration 
of chiral symmetry in the chiral limit. 



> 





Fig. 1. Angle variation of for different values of tempera- 
tures and chemical potentials in the case of chiral limit. The solid 
line corresponds to T = 150MeV,yU = lOOMeV, dashed line cor- 
responds to T = 250MeV,yU = lOOMeV, dash-dotted line cor- 
responds to T = 40MeV,yU = 300MeV, and dot corresponds to 
T = 150MeV !y u = 300MeV. 
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Fig. 3. The conventional chiral condensate -{cr) n and the dressed 
Polyakov loop 27] as functions of temperature for different values 
of the chemical potentials in the chiral limit. Here, -{cr) n and E\ 
both are measured in [GeV 3 ]. From right to left the values of the 
chemical potential are 0, 200, 300MeV, respectively. 



Almost no variation with respect to angle is found for 
low temperatures. As the temperature increases the varia- 
tion over the angle grows. We expect the absolute value of 
the chiral condensate decreases with the increase of tem- 
perature. However, from the figure, this conventional be- 
havior of the chiral condensate with temperature only per- 
sists up to a certain angle, beyond which the opposite be- 
havior is observed. The plateau around <p = n is more flat 
above T c in case of zero current quark mass and its value 



Fig.[3]and|4]show the behavior of the conventional chi- 
ral condensate -{cr) n and the dressed Polyakov loop S\ 
at different chemical potentials as functions of tempera- 
ture for m — and m = 5.5MeV, respectively. For both 
cases, it is observed there are three temperature regions 
for —{<t)„ and E\. For —{&)„, at smaller temperatures it 
remains constant at a value corresponding to the value of 
the conventional chiral condensate in the vacuum, then it 
rapidly decreases in a small window of temperature and 
eventually almost saturates to a lower value. The rapid de- 
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Fig. 4. The conventional chiral condensate -{cr) n and the dressed 
Polyakov loop £ { as functions of temperature for different val- 
ues of the chemical potentials in the case of explicit chiral sym- 
metry breaking m = 5.5MeV and —{cr) n , 2^ are measured in 
[GeV 3 ]. From right to left the values of the chemical potential 
are 0, 200, 300MeV, respectively. 
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Fig. 5. Two-flavor phase diagram in the T - yu plane for the case 
of chiral limit. The solid line is the critical line for Z\, and the 
dashed line is the critical line for conventional chiral phase tran- 
sition. The solid circle indicates the critical point for chiral phase 
transition. 



creasing occurs at different temperatures for different val- 
ues of the chemical potentials. On the other hand the be- 
havior for the dressed Polyakov loop is just the opposite. It 
remains almost zero for small temperatures and then rises 
rapidly, finally saturates to a high value which varies very 
slowly with temperatures. The almost zero value of E\ for 
small temperatures is due to the fact that the U(\) bound- 
ary condition dependent general quark condensate nearly 
does not vary with the angle (f> for small temperatures (see 
Eq.©. 

For finite quark mass, near the critical temperature re- 
gion, both -(cr) ff and E\ change more slowly than those in 
the case of chiral limit. 

Fig. [5] and [6] show the T - /v phase diagram for the 
case of m = and m = 5.5MeV, respectively. The tran- 
sition temperatures are calculated from the slope analysis 
of the conventional chiral condensate {<r) n and the dressed 
Polyakov loop. The transition temperatures calculated from 
the conventional chiral condensate represent the chiral phase 
transition temperature. On the other hand the behavior of 
the dressed Polyakov loop is supposed to indicate the de- 
confinement transition temperature. In our present frame- 
work confinement is not accounted for. However, if we 
look at the curves presented in figure [3] and |H they still 
show an order parameter like behavior. We would like to 
point out here that for the same reason as stated above 
we are not concerned here with the order of the phase 
transtion from the dressed Polyakov loop and all the com- 
ments about the order of the phase transition below are 
with respect to the chiral phase transition only. 

For m = case, we find almost exact matching for the 
transition temperatures calculated from these two quanti- 
ties in the whole T - ft plane as shown in Fig. [5] 

For the case of finite quark mass m - 5.5MeV, it is 
observed from Fig.|6]that the two critical temperatures are 



0.24 




0.1 0.2 0.3 

H [GeV] 

Fig. 6. Two-flavor phase diagram in the T - \i plane for the case 
of m = 5.5MeV. The solid line is the critical line for S\, and the 
dashed line is the critical line for conventional chiral phase tran- 
sition. The solid circle indicates the critical end point for chiral 
phase transition. 



different in the low baryon density region. The difference 
however decreases from low to high chemical potential, 
and the two critical temperatures start to match around 
the critical end point for chiral phase transition. For zero 
chemical potential and small current quark mass m = 5.5MeV, 
we find about 7MeV difference between T x c and T®, and 
T* < T®. Similar trend has been observed in another study 
based on Dyson-Schwinger approach P31l . where they found 
chiral transition to occur about 10 - 20MeV below the 
deconfinement transition. Though these studies are not a 
complete one and these differences may be due to the ef- 
fects of crossover transition. As pointed out in l25ll during 
crossover, different observables are expected to behave dif- 
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ferently and there is no way to define a unique crossover 
temperature. 

We extend our study further to see what happens if we 
increase the current quark mass further. We find at zero 
chemical potential, the difference between the transition 
temperatures calculated from dressed Polyakov loop and 
conventional chiral condensate increases as we increase 
the current quark mass (see Fig. |7j. Initially the difference 
is zero for zero current quark mass but for m = 200MeV 
we find about 26MeV difference between the two temper- 
atures. It is worthy of mentioning that this result is just 
for illustrative purpose as there are limitation of using NJL 
model with such a huge current quark mass. 
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Fig. 8. Temperature derivative of the general chiral condensate 
-jjr- for m = and li = at three temperature cases: T < T c x , 
T = T c x and T > T c x , where T c x is the chiral transition tempera- 
ture. 



finite current quark mass) and the region with large values 
of d(cr)j,/dT shrinks. Therefore, the integral in Eq. (0 i.e. 
dZildT decreases as temperature increases. In all, dE\jdT 
gets its maximum at T x and the two transition tempera- 
tures T® and T x coincide in the case of second order chiral 
phase transition. 



Current quark mass [MeV] 

Fig. 7. The critical temperatures T x and T® for different current 



4 The relation between T x c and T® 

In the following, we offer a possible understanding on the 
simultaneity of the transition temperatures for 1st and 2nd 
order chiral phase transitions and the apparent difference 
between the two for the case of crossover. 

As mentioned earlier the transition temperatures are 
determined from the slope analysis of the conventional chi- 
ral condensate and the dressed Polyakov loop. So let us 
look at the temperature derivative of the general chiral con- 
densate d{cr)^ldT as functions of <p, the integral on which 
gives the temperature derivative of the Dressed Polyakov 
loop 

dT 

Fig. [8] shows d{cr)^ldT at different temperatures in the 
cases of second order chiral phase transitions. Around T x , 
large values of d{cr)^ldT appear around <p = n and domi- 
nate the integral in Eq. (0. Below T x , d(cr)$/dT around 
<f> — n increases monotonously as temperature increases. 
As a result, dZ\ldT increases. Above T x , d{cr)$ldT in the 
center region becomes zero (or very small in the case with 



dz { _ r 2n d(t> 

dT Jo 2;t 



(7) 




Fig. 9. Temperature derivative of the general chiral condensate 



for m 



dT 

T < T x , T = 
temperature. 



= and /J = 300MeV at three temperature cases: 
T c x and T > T c x , where T c x is the chiral transition 



In the case of first order chiral phase transition, the situ- 
ation is more complicated. Due to the discontinuity of (cr)^ 
at the first order chiral phase transition point, the temper- 
ature derivative of the dressed Polyakov loop can be ex- 
pressed as 



dZ\ _ r 2 "d(P 
dT Jo 27r 



dT 



Cos0 c d(f> c 

d{o-) c — , 

7T dl 



(8) 
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where, the first term is determined by the regular behavior 
of d(cr)$/dT (see Fig|9]l, the second term is due to A{cr) c , 
the jump of (cr)^ at the first order phase transition point at 
<p = (p c . When T < T x , the second term vanishes. Now 
we consider two limiting cases. First, in the case of a weak 
first order phase transition, A{cr) c is small and d(cr)^/dT 
around cp — n is large, as showed in Fig(9] So the first term 
dominates the result of Eq. © and gives the similar result 
as that in the case of a second order chiral phase transition. 
Second, in the case of a strong first order chiral phase tran- 
sition, A{cr) c is large and d{cr)^ldT is small. So the second 
term dominates the result of Eq. ©. Then is strongly 
dependent on the detailed information of d(p c /dT. Our nu- 
merical results show that d(p c /dT decreases as temperature 
increases, and so the second term also gives a decreasing 
contribution. In all, it is clear that dX[/dT gets its maxi- 
mum at T x , i.e. T® and T x coincide in the case of a weak 
first order chiral phase transition due to remnants of second 
order chiral phase transition. The coincidence in the case 
of a general first order chiral phase transition is supported 
by our numerical results and can be generally expected. 




-0.3 1 1 1 1 1 

0.5 1 1.5 2 

(|)jt 

Fig. 10. Second derivative of the general chiral condensate -j^r 
for m = 5.5MeV and \i = at three temperature cases: T < 
T c x , T = T c x and T > T c x , where T/ is the chiral transition 
temperature. 

For the case with finite quark mass and small chemi- 
cal potential, we have no phase transitions but crossover. 
So let us consider the second temperature derivative of the 
dressed Polyakov loop 

dT 2 Jo 2tt dT 2 ' W 

whose zero point determines the transition temperature T®. 
Fig. [10] shows the second derivatives of the general chiral 
condensate d 2 (cr)^/dT 2 at three temperature cases: T < 
T c x , T = T c x and T > T/, where T c x is the chiral transi- 
tion temperature. Similar to our previous observation, large 
values of d 2 (cr) ( f,/dT 2 appear around <p — it (see the two 
maximums in Fig. [Tol l, as remnants of the second order 



phase transition in chiral limit. The difference is that d 2 (o-),f,/dT 2 
in the center region (see the minimum in Fig. QJJi, is sup- 
pressed below T c x , approaches zero at T c x and changes its 
sign above T c x , For T < T c x , d 2 {(r)$ldT 2 around the two 
maximums dominate the integral in Eq.© and d 2 E\ldT 2 
does not change its sign. Above T c x , the negative part around 
the minimum cancels the contributions from the maximums, 
and up to a certain temperature T®, this cancelation leads 
to the zero of the integral in Eq.®- In all, the zero point 
of d 2 Z\ldT 2 comes from the negative contribution of the 
minimum at T > T c x , so the transition temperature related 
with the dressed Polyakov loop must be higher than the 
chiral transition temperature, i.e. T L P > T c x . 



5 Phase diagram for three-flavor case 

We now consider the phase diagram for three-flavor case, 
and the parameters are taken from Ref. 115 11521 : 





m„[MeV] 


m s [MeV] 


G s A l 




Chiral-limit 








1.926 


12.36 


finite-mass 


5.5 


140.7 


1.918 


12.36 



Table 1. Two sets of parameters in 3-flavor NJL model: the cur- 
rent quark mass m q for up and down quark and m s for strange 
quark, coupling constants G and K, with a spatial momentum cut- 
off A = 602.3 MeV. 



In the chiral limit, i.e, for the case of m u — nid — m s — 
case, we find almost exact matching for the transition 
temperatures calculated from these two quantities in the 
whole T - [x plane as shown in Fig. Q~TJ The chiral phase 
transition is of 1st order in the whole T - y, plane. 

For the case of finite quark mass m u = = 5.5MeV 
and m, = 140.7MeV, it is observed from Fig. Q~2]that the 
the chiral and deconfinement phase transitions are flavor 
dependent. For the degenerate light u, d quarks, the result 
of phase transitions are similar to that in Fig. [6] The critical 
temperatures of chiral and deconfinement phase transitions 
have very small difference in the low baryon density re- 
gion when the transition is of crossover, and the difference 
vanishes at the critical end point. Here the CEP is located 
at (T^ P ,fj^ p ) = (91MeV,315MeV), which is different 
from Fig. [6] This difference comes from: 1) different model 
parameters have been used, 2) the coupling of s quark to 
u, d quark contributes one extra term in the thermodynam- 
ical potential comparing with the pure two-flavor case. For 
the relative heavy s quark, it is found that this quark has 
separate phase transition lines. At low baryon density, the 
variation of s quark condensate in the crossover region is 
so diffused that it is not possible to identify the crossover 
temperature. When density increases, the crossover tem- 
peratures for chiral and deconfinement can be effectively 
extracted and still shows the relation of T x < T®, These 
two critical temperatures of s quark coincide at the critical 
end point (T s CEP ,n s CEP ) = (75MeV,445MeV). 
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Fig. 11. Three-flavor phase diagram in the T — ft plane for the 
case of chiral limit. The solid line is the critical line for £ t , and 
the dashed line is the critical line for conventional chiral phase 
transition. 




Fig. 12. Three-flavor phase diagram in the T-fi plane for the case 
of m u = nid = 5MeV and m s = 140.7MeV. The solid lines are 
the critical line for L\ , and the dashed lines are the critical line for 
conventional chiral phase transition. The solid circles indicate the 
critical end points for chiral phase transitions of u, d and s quark, 
respectively. 



6 Conclusion and discussion 

We investigate the chiral condensate and the dressed Polyakov 
loop or dual chiral condensate at finite temperature and 
density in the two and three-flavor Nambu-Jona-Lasinio 
model. We find the behavior of dressed Polyakov loop in 
absence of any confinement mechanism still shows an or- 
der parameter like behavior. It is found that in the chiral 
limit, the critical temperature for chiral phase transition 
coincides with that of the dressed Polyakov loop. In the 
case of explicit chiral symmetry breaking, it is found that 
the phase transitions are flavor dependent, and the critical 
temperature for chiral transition T x is smaller than that of 
the dressed Polyakov loop T® in the low baryon density re- 



gion where the transition is a crossover. With the increase 
of current quark mass the difference between the two crit- 
ical temperatures is found to be increasing. However, the 
two critical temperatures coincide in the high baryon den- 
sity region where the phase transition is of first order. For 
three-flavor case, it is observed that there are two critical 
end points in the (T,/j) phase diagram. 

From symmetry analysis, the dressed Polyakov loop 
can be regarded as an equivalent order parameter of de- 
confinement phase transition for confining theory. In the 
NJL model, the gluon dynamics is encoded in a static cou- 
pling constant for four point contact interaction. Since in 
this work we have included only quark degrees of free- 
dom, a quantitative comparison will not match with other 
results. But the interesting fact is that the qualitative fea- 
tures (angle variation, temperature variation) of the dressed 
Polyakov loop remains the same. Moreover, we expect that 
independent of the input of gluedynamics to the quark prop- 
agator, it is a general feature for T x = T® in the case of 1st 
and 2nd order phase transitions, and T x < T® in the case 
of crossover, which qualitatively agrees with the lattice re- 
sult in Ref. [25 1. This might indicate that for full QCD, in 
the crossover case, there exists a small region where chiral 
symmetry is restored but the color degrees of freedom are 
still confined. This result should be checked in other ef- 
fective models, e.g. in the framework of Dyson-Schwinger 
equations (DSE). 

The (T,/j.) phase diagram with three flavors and with 
Ua(1) anomaly as well as diquark condensate will be stud- 
ied in the near future. 
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